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Introduction 

Knowledge of the oxidation state of a magma is critical as it is one of the parameters which controls the 
nature and composition of the resulting crystals [Eugster, 1957; Eugster and Wones, 1962; Carmichael, 
1967; Frost and Lindsley, 1991]. In terrestrial magmatic systems, oxygen fugacity (f 02 ) is known to vary 
by over nine orders of magnitude [Carmichael, 1991; Parkinson and Arculus, 1999]. With variations of this 
magnitude, understanding the compositional differences, phase changes, and crystallization sequence 
variations, caused by the magma f 02i is essential in deciphering the origin of all igneous rocks. 

Magmatic oxidation state is of great importance in that it reflects the degree of oxidation of the source 
region and can provide insight into magmatic processes, such as metasomatism, degassing, and 
assimilation, which may have changed them. Carmichael [1991] argues that most magmas are unlikely to 
have their redox states altered from those of their source region. This assumption allows for estimation of 
the oxidation state of planetary interiors. Conversely, it is known that the / 02 of the magma can be affected 
by other processes, which occur outside of the source region and therefore, the oxidation state may record 
those too. Processes which could overprint source region fugacities include melt dehydrogenation or other 
volatile loss, water or melt infiltration, or assimilation of oxidized or reduced wallrock [Brandon and 
Draper, 1996; Ague, 1998; Frost and Ballhaus, 1998; Parkinson and Arculus, 1999]. Understanding which 
of these processes is responsible for the redox state of a magma can provide crucial information regarding 
igneous processes and other forces active in the region. 

The composition of the SNC basalts and their widely varying proposed oxidation states raise some 
interesting questions. Do the SNC meteorites have an oxidized or reduced signature? What was the 
oxygen fugacity of the SNC source region at the time of melt generation? Is the fugacity calculated for the 
various SNC samples the fugacity of the magma source region or was it overprinted by later events? Are 
there different oxidation states in the Martian interior or a single one? This proposal seeks to address all of 
these questions. 

Progress Report 

Task la of the original proposal has been completed. All present publications resulting from this are 
presented in the references. This task consisted of an assessment of the variation of Eu- valence state in 
early crystallized, experimental low-Ca pyroxenes (pigeonite, Pgt) with magmatic f 0 2 , with an emphasis on 
constraint of the oxidation state recorded in the basaltic shergottites. Fugacity was measured in the 
pyroxenes using two techniques which utilize oxidation state proxies of D(Eu/Gd) pgt / me j 1 and 
D(Eu/Sm) pgt/me i t . The results of this task provide a new way to measure Martian f 02 . 

Task 1 . SNC fo 2 recorded in pyroxenes: An experimental study. 

Experiments to characterize pyroxene composition have been conducted over a range of 
representative magmatic f 02 s. Dry crystallization experiments were performed on a starting composition 
composed of a synthetic shergottite melt in equilibrium with liquidus pyroxenes, A* M (see Johnson et al. 
[1991]). The melt was doped with ~lwt.% each of Eu 2 0 3 , Gd 2 0 3 , and Sm 2 0 3 following the method of 
McKay et al. [1986]. Experiments were run at the Johnson Space Center in one-atmosphere, C0-C0 2 gas 
mixing furnaces using Re-loops. The dry experiments were run over a range of representative magmatic 
fo 2 s 9 from highly reduced IW-2 to moderately oxidized IW+3.5 (equivalent to QFM, quartz-fayalite- 
magnetite). The f Q2 was controlled through the mixing of the CO and C0 2 gases, and measured with an 
external Zr0 2 sensor. All dry runs were held at 1235°C for 30 minutes to ensure melt homogenization. 

Runs were then cooled to 1 165°C, 15-20°C below the dry liquidus temperature determined by Minitti and 
Rutherford [2000] for a similar starting composition, at 60°C/hour. Charges were held for 24 hours to 
facilitate crystal nucleation and growth. Samples were then rapidly quenched and the products thin 
sectioned. 

Hydrous crystallization experiments were also performed on A* M in TZM (tungsten-zirconium- 
molybdenum) pressure vessels at Brown University. Experiments were run at 1 1 10°C and Ptotar'Pmo^OO 
bars. Experimental f Q2 s ranged from oxidizing, nickel-nickel oxide (NNO = IW+4.5), to reducing, 
graphite-methane (fW+0.5). The higher f 02 runs were buffered using pre-determined mixtures of Ar and 
CH 4 gases appropriate for the desired f Q2 . These runs included powdered buffer assemblages in separate 



tubes within the experimental capsules to monitor the f 02 . For example, Ni (m ) and NiO were placed inside 
all NNO runs. At the end of the run the presence of both phases was verified. Runs which contained only 
one phase were discarded. The low f Q1 experiment was buffered through the breakdown of CH 4 gas. To 
better approximate the liquidus, runs were initially pressurized to 500 bars and held for 24 hours. Runs 
were then slowly decompressed over ~2 hours to a final pressure of 200 bars, corresponding to 
approximately 1.5 wt.% water, an amount determined to be plausible for SNC magmas by Minitti and 
Rutherford [2000]. This pressure drop initiated pyroxene crystallization. Runs were held for at least 24 
hours to facilitate crystal nucleation and growth. Samples were then rapidly quenched and the products 
thin sectioned. 

All experimental phenocryst phases and melts were analyzed for their major element compositions on 
the Brown University Cameca Camebax and SX-100 electron microprobes. Eu/Gd and Eu/Sm ratios in the 
pyroxenes and coexisting melt were measured on the Woods Hole Oceanographic Institution Cameca 3F 
and Arizona State University Cameca 6F ion microprobes. 

The results of both the dry and hydrous experiments are given in Figure 1 [McCanta and Rutherford, 
2001; McCanta et al., 2002a]. For a dry basaltic shergottite melt, the average D(Eu/Gd) pyx/me i t range from 
0.156(/o2 = rW-l) to 0.639(QFM). D(Eu/Sm) pyx / me jt range from 0.279(IW-1) to 1.135(QFM). The 
D(Gd/Sm) pyx / me it is 1.76±0.02, indicating all variation in Figure 1 is due to Eu variation. This new Eu- 
oxybarometer calibration for pigeonite yields refined / 02 values for the basaltic shergottites, ranging from 
IW+3.1 for Shergotty to IW-fO.2 for QUE94201. These numbers differ in value, but not trend, from those 
published previously [e.g. Herd et al., 2001; Wadhwa, 2001]. These f 02 values were calculated using 
natural pigeonite core data from Wadhwa et al. [1994]. The f 02 s calculated from our pigeonites [McCanta 
et al., 2002b; McCanta and Rutherford, 2002] are similar to those determined in a recent calibration of 
augite REE distribution coefficients vs. f 02 [Musselwhite and Jones, 2002; Musselwhite and Wadhwa, 
2002] and to those determined using recent Fe-Ti oxide data [Herd et al., 2000] (Figure 2A and B). 

For a hydrous basaltic shergottite melt the average D(Eu/Gd) pyx / melt range from 0.62 (/^IW+O^) to 
0.69 (IW+4.5). D(Eu/Sm) pyx / melt was not measured in hydrous runs as the starting material was only doped 
with Gd 2 0 3 and Eu 2 0 3 . The higher f 02 results determined in the dry experiments correspond well with 
results from the hydrous experiments. 

Our experimental data provides a calibration of the effects of f 02 on pigeonite REE partitioning in both 
dry and hydrous basaltic shergottite melts. Because pigeonite is one of the first phases to crystallize, this 
technique allows for the measurement oifoi during the earliest crystallization history of the melt. Both 
distribution coefficients (Eu/Gd and Eu/Sm) record comparable f 02 values for measured Martian meteorites 
indicating that either method can be used in future analyses. 

Our data suggests that the Martian meteorites crystallized under a variety of different f 02 conditions. 
This has been recognized previously in studies of augite [e.g. Wadhwa, 2001]. Our new results from the 
pigeonite Eu-oxybarometer detailed above show ~2-3 orders of magnitude difference in f Q2 between the 
basaltic shergottite meteorites, not as large a range as previously argued. This variation is less than the 
observed terrestrial variation of ~9 orders of magnitude [e.g. Carmichael, 1991], but much greater than that 
seen on the moon. Several possibilities exist that could explain the range in Martian f 02 . It could reflect 
intrinsic f 02 differences in the magma source regions on Mars (i.e. the Martian mantle or lithosphere), or it 
could result from interactions between the magma and an assimilant, that is as yet unsampled, either at 
depth or near the surface. This new f 02 data may provide insight into which of these is the more likely 
explanation. 

Experimental studies have shown that pigeonite is one of the first minerals to crystallize from a 
shergottite melt, whereas Fe-Ti oxides come on the liquidus late in the crystallization sequence (after 30 
and 60% crystallization in hydrous and dry experiments, respectively, to possibly 95% in Los Angeles) 
[e.g. Minitti and Rutherford, 2000; Xirouchakis et al., 2002]. Therefore, pigeonite core f 02 values should 
record a different portion of magmatic history than the oxides. However, when the pigeonite f 0 2 values are 
compared with the most recent determination of Fe-Ti oxide f 02 values, a close correlation is observed 
(Figure 2B). This implies that the melt f 02 for each of these meteorites was buffered internally throughout 
much of the crystallization sequence. That is, the early crystallizing minerals were exposed to the same 
relative f 02 environment as the late crystallizing minerals. These results are consistent with the argument 
that the variations in f 0 2 recorded in the basaltic shergottites are the result of intrinsic f 02 differences in the 
magma source regions of these rocks. However, assimilation cannot be completely ruled out on this basis 



alone. If the assimilation occurred before any crystallization had taken place, the primary ^02 recorded in 
the phenocrysts would be that of the source region modified by the assimilant. 

Task lb: Investigation of Henry f s Law behavior of REE in pigeonite. 

Several tasks were originally proposed to help illuminate the oxidation state of the Martian interior. 

The first task, discussed in the progress report above, has been completed as stipulated in the original 
proposal. An additional question arose as a result of the data collected in Task la relating to the Henry’s 
Law behavior of the REEs which were doped to weight percent levels in the experiments. In order to 
address this question, new experiments were completed on starting compositions containing 0.5 and 0. 1 
weight percent REE, respectively, as opposed to the original starting composition doped with 1 .0 weight 
percent REEs. These runs were run in the gas mixing furnaces at Johnson Space Center over the same 
range of /02 conditions as the original experiments. All run products were analyzed for their REE 
concentrations on the Arizona State University Cameca 6F ion microprobe. The data collected was 
compared to the data presented above and it was determined that Henry’s Law had not been violated in any 
of the Eu-oxybarometer calibration experiments (Figure 3). This task provided confidence in the pigeonite 
Eu-oxybarometer and the Martian crystallization f 02 values the pigeonite record. 
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Figure Captions 

Figure 1. Dry and hydrous calibration of D(Eu/Gd) pigeonite versus f 02 (in log units relative to the IW buffer). 
Error bars on the distribution coefficients indicate the standard deviation between ion microprobe analyses 
of the materials in the same sample. Solid line is the best- fit curve (2 nd order polynomial) for each system. 

Figure 2. (A) Comparison with augite Eu-oxybarometer calibration of Musselwhite and Jones (2002) 
using shergottitic melt starting composition. S: Shergotty, EA: EETA79001 lithology A, EB: EETA79001 
lithology B, Z: Zagami, Q: QUE94201. Pigeonite f 02 error bars are 2a. Augite fo 2 error from Musselwhite 
and Jones (2002). (B) Comparison of f 02 calculated using dry pigeonite Eu-oxybarometer versus Fe-Ti 
oxide equilibria (Herd et aL> 2001). Solid line is a 1 : 1 line. Pigeonite f 02 error bars are 2a. Oxide f 02 error 
from Herd et al (2001). 

Figure 3. Dry calibration of D(Eu/Gd) p j ge omte/meit and D(Eu/Sm)p igconitc/m eit versus f 02 (in log units relative to 
the IW buffer) for 1.0 wt.% REE experiments. Pgt: pigeonite. Solid line is the best-fit curve (2 nd order 
polynomial) for each system. Squares are measured data for Henry’s Law experiments with 0.1 wt.% of 
REEs. Note that the low-REE runs plot either on top of or within error of the high-REE runs. For 
example, at f 02 = IW+3.5, the two D(Eu/Gd)pi gconit e/ ra eH data points are identical, while the two 
D(Eu/Sm) p igeonjte/meit data points are the same within error. Error bars on the distribution coefficients 
indicate the standard deviation between ion microprobe analyses of the materials in the same sample. 
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